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The region of small transverse momentum in qq— and gg— initiated processes must be studied 
in the framework of resummation to account for the large, logarithmically-enhanced contributions 
to physical observables. In this letter, we study resummed differential cross-sections for Higgs pro- 
duction via bottom-quark fusion. We find that the differential distribution peaks at approximately 
15 GeV, a number of great experimental importance to measuring this production channel. 
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Resummation of total and differential cross-section for 
the inclusive production of a Higgs boson has concen- 
trated on t he g luon-gluon initial stated IS IS IS IS IS 13 
IS IS 13 IH Ea USUI- In the Standard Model (sm), the 
gluon-gluon initial state gives the largest contribution to 
the total and differential cross-sections, but this is not al- 
ways the case in extensions of the SM. In the Minimal Su- 
persymmetric Standard Model (mssm) the bottom-quark 
fusion initial state can be greatly enhanced, perhaps lead- 
ing to the observation of a supersymmetric signal in na- 
ture, if the location of the peak of the differential distri- 
bution in known. 

The MSSM contains two Higgs doublets, one giv- 
ing mass to up-type quarks and the other to down- 
type quarks. The associated vacuum expectation val- 
ues (vevs) are labeled v u and va respectively, and fix 
the free MSSM parameter tan/3 = v u /vd- In the MSSM, 
there are five physical Higgs boson mass eigenstates. In 
this letter, we are interested in the neutral Higgs bosons 
{h°, H°, A } which we will call <E> generically. 

In contrast to the SM, the bottom-quark Yukawa cou- 
plings in the MSSM can be enhanced with respect to 
the top-quark Yukawa coupling. In the SM, the ratio 
of the ti& and 66$ couplings is given at tree-level by 
Aj M /A| M = m t /m b ps 35. In the MSSM, the coupling de- 
pends on the value of tan (3. At leading order, 
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where a is the mixing angle between the weak and the 
mass eigenstates of the neutral scalars. Given the mass 
of the pseudoscalar M^o and tan j3, the angle a can be 
determined given reasonable assumptions for the masses 
of the other supersymmetric particles in the spectrum. 

The form of /$ shows us that the production of the 
pseudoscalar due to bottom-quark fusion is enhanced by 
a factor of tan 2 /?, which is a free parameter in the theory. 



THE BOTTOM-QUARK 

It is also important to define what is meant by a 
bottom-quark distribution^^ ^S El- In our analy- 
sis, we employed the CTEQ6.1M bottom-quark parton 
distribution | lSj with a s (Mz) — 0.118 and set the mass 
of the Higgs boson of interest M$ = 120 GeV. We com- 
pared the bottom-quark distribution function in the PDF 
set and the numeric solution to the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (dglap) equations for a single 
quark splitting. A bottom-quark distribution for a gluon 
splitting into a bb pair can be written in the DGLAP for- 
malism as 



g(y,v), (3) 



where g is the gluon distribution, and the DGLAP splitting 
function is 



P q ^ g (z) = l[z 2 + (l-z) 2 }. 



(4) 



The bottom-quark distribution is encoded into the 
CTEQ PDF set in this manner^^ , but takes into account 
multiple quark splitting functions. When evaluated with 
the dglap formalism, we found the differential cross- 
section increased by approximately 10% near its peak. 
However, we used the native bottom-quark distributions 
for speed and to understand their built-in uncertainties. 

Previouslv[l^. we calculated in detail the resumma- 
tion coefficients for a differential cross-section for the 
scalar and pseudoscalar Higgs boson from the gluon- 
gluon initial state. In this letter, we will calculate the 
resummation coefficients needed for the resummation of 
the bb initial state for the scalar and pseudoscalar Higgs 
bosons in the same manner as the gluon-gluon channel in 
Ref. 0| ■ We will leave the bottom-quark-Higgs coupling 
set equal to the SM value so that the reader can scale the 
results to whatever coupling value is of interest. 
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RESUMMATION 

The resummation formalism needs the lowest order to- 
tal cross-section as a normalization factor (see fl4| for 
details), bb — > <I> in this case. Following Ref. |20j, we will 
ignore the bottom-quark mass except in the Yukawa cou- 
pling with the Higgs boson. Although the pseudoscalar 
Higgs couples to quarks with a 75, there are no differ- 
ences in the matrix elements (modulo the MSSM coupling 
factor) when the bottom-quark mass is neglected. 

It is important to use the MS running mass for the 
bottom-quark in our calculation as the difference from 
the pole mass at the scales involved is considerable |2ll 
122^ . In the SM, the bottom-quark Yukawa coupling is 
\f* = sflSm^Jv, where v is the SM vev and is approxi- 
mately equal to 246 GeV and m& is the MS running mass. 
We have set the bottom-quark mass TOb(m&) = 4.62 GeV 
in our calculations. The NLO running of the bottom- 
quark mass corresponds to mfe(M$) = 3.23 GeV. The 
coupling in the MSSM can be written 
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The spin- and color-averaged total partonic cross- 
section (see Fig. QJi) for the leading order subprocess, 
bb — > $, can be easily written 
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where z — M^/s and the number of colors N c = 3. 
We also need the LO differential cross-section (Fig. QJ:) 
for the next-to-leading log (NLL) resummation coeffi- 
cients for the differential cross-section. If we remove the 
<5(1 — z) factor from our prefactor cto, then we can write 
the spin- and color-averaged differential cross-section for 
Kpi)Kp2) 9(-P3)$(-P5) as 
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where C F = (iV 2 — 1)/2N C , and the kinematic variables 



are defined as s = (P1+P2) , t = (pi+Ps) 



(P2+P5) , 



and M| = p£- In our second line, we have written the 
differential cross-section in terms of ut = p^s for the 2 — > 
2 process. 



To find the resummation coefficients for a differential 
cross-sections (3, 0, 0, El we integrate the differential 
cross-section around p t = 
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and label this result 'real' as it is similar to the real 
corrections to the LO total cross-section. Working in 
N = 4 — 2e dimensions, we find 
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To regularize this result, we need to add the virtual 
corrections that are shown in Fig. ^3. These corrections 
are very similar to Drell-Yan corrections [2^. The virtual 
corrections can be written as 
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-C F factor would be 



-AC 



F m 



In the Drell-Yan case, the 
When the two results are added together the resumma- 
tion coefficients are easily read off from the expression. 
The total expression is 
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Keeping with the notation of Ref. [I 
coefficients with an overbar as follows 



we write these 
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It is important to note that in contrast to W^/Z® pro- 
duction and Drell-Yan processes [13, [3], the C^ 1 ** coeffi- 
cient is positive. 

Finally, let us turn to determining the C^ and C' 2 ' 
coefficients for the total cross-section resummation, al- 
though total cross-sections will not be presented in this 
letter. Using the results of Ref. [iof. we take the Merlin 
moments of the corrections in the limit N — > 00 (z — > 1 
in z-space). The NLO corrections are easy to color de- 
compose due to the presence of only one color factor [Til]. 
Leaving the terms that were originally proportional to 
the 5(1 — z) factor inside curly brackets, we find 
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+ |2C F C2-C F + 21n^-|+2C F C2 (13) 

where we have given both the renormalization scale jj, r 
and the factorization scale [if dependence in the results. 

In contrast to NLO, the NNLO corrections contain a 
mix of color factors (both C A and C F appear) . Although 
it is easy to see that the factor proportional to ln 4 (iV) 
should clearly be 2C F , no unique color decomposition 
from the results provided in Ref. |2lj can be determined 
for all the terms in the expression. However, the numeric 
result can be written 




(a) (b) (c) 

FIG. 1: Diagrams needed for the bb initial state resummed dif- 
ferential cross-section. Figure^, is the lowest order produc- 
tion channel and couples differently for different Higgs bosons. 
Figure^} is the virtual correction to the lowest order process. 
Figure 0; is the lowest order graph contributing to the differ- 
ential cross-section. The crossed graph is not shown. 
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We can also determine the NNLL A^ and -B^ 2 ) coef- 
ficients, agrees with a previous calculation}^. We 
find 
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for the -E^ 1 ) and B^ coefficients. The Mcllin moments 
A^' and are novel, as is B?'. 

bb bb b 



RESULTS AND CONCLUSIONS 

The differential resummation coefficients and the po- 
sition of the peak of the differential cross-section is of 
great interest to the experimental community involved 
with Higgs research at the LHC, particularly in the 
71/$ = 120 GcV mass range. Here the Higgs will de- 
cay primarily into 66 pairs that can be tagged. Knowing 
where the peak of the differential distribution lies, espe- 
cially if it is below the pt of a typical trigger event, is of 
utmost importance. This letter will help in the analysis 
of the 66 initial state. 

The results of our calculations can be found in Fig- 
ure^ We have done our analysis for the LHC (a proton- 
proton collider at VS = 14 TcV). Wc find that the dif- 
ferential distribution at the LHC peaks at a transverse 
momentum of approximately 15 GcV. We find that the 
magnitiude of the differential cross-section is an excel- 
lent match with previously published results |2ll |25| . 
The results for the Tevatron are extremely similar, but 
are smaller by a factor of 60 and the peak moves to a 
transverse momentum of approximately 13 GeV in the 
differential distribution. 

A detailed study of the uncertainties in the calculation 
show that the uncertainty due to the PDF set is approx- 
imately 8 — 12%. At the peak of the distribution, the un- 
certainty is approximately 10% due to the PDFs. When 
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FIG. 2: Figure |2t shows the errors associated with the 
CTEQ6.1M PDF set. The variation is approximately 8-12%. 
Figure^ show the variation of the renormalization and fac- 
torization scale for a factor of 1/4 and 4. These scales were 
chosen because there has been great interest in the scale 
fi = M$/4. We find this variation to be approximately 20%. 
Figure [5Jj shows the movement of the peak of the differen- 
tial distribution (at 15 GeV) for a variation in the scale by a 
factor of 10. 



the scale is varied by a factor of four, we see a variation in 
the differential cross-section of approximately 20%. This 
would give us a combined uncertainty of 32%, which is 
slightly better than the gluon-gluon channel[14J uncer- 
tainty in the differential distribution. However, when the 
scale is only varied by a factor of two (as was the case for 
the gluon-gluon channel), the total uncertainty drops to 
25%. 

We have calculated the rcsummation coefficients 
needed for NLL inclusive Higgs production via bottom- 
quark fusion in the SM and the MSSM for the differential 
cross-section and for the NNLL resummation for the to- 
tal cross-section. We find a smaller uncertainty in the 
bottom-quark initial state than the gluon-gluon initial 
state. 
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gelsong, F. Olness, and A. Ficld-Pollatou. 1 would also 
like to thank W. Kilgore and R. Harlander for supply- 
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